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Abstract: We continue to investigate possible signatures of a pre-inflationary anisotropic 
j — phase in two-point and three point correlation functions of the curvature perturbation for 

high-momentum modes which exit the horizon well after isotropization. The late time dy- 
namics of these modes is characterized by a non-Bunch Davies vacuum state which encodes 
^vq all the information about initial anisotropy in the background space-time. We observe that, 

unlike the non-planar momenta, there exist regimes of planar momenta for which scale in- 
variance of the power spectrum is strongly broken. This regime of planar momenta gives rise 
^ to enhanced non-Gaussianity in certain squeezed triangle configurations, although the en- 

hancement of the Jnl parameter is limited by the breakdown of linear perturbation theory at 
"exact planarity" . Finally, we demonstrate that for the range of planar modes for which scale 
invariance of the power spectrum is preserved, non-Gaussianity in the curvature perturbation 
spectrum is naturally constrained to be extremely small. 
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1. Introduction 

In addition to the power spectrum, understanding detailed features of the inflationary sce- 
nario requires probing three-point correlation functions and beyond, collectively referred to as 
non-Gaussianitiy [1]. It is well-known that observable non-Gaussianity requires a departure 
from the standard single-field inflation with a canonical action [2, 3]. Multiple scalar fields, 
non-canonical action for the scalar field, introducing higher derivative terms in the action or 
having a non-standard vacuum state (see [4] for reviews) are some of the variants from the 
standard inflationary model that have been studied. A common feature of all these models, 
however, is that they have a homogeneous and isotropic background for the perturbations to 
evolve. 

Investigating possible signatures of a strongly anisotropic phase of expansion at early times 
in the late time cosmological perturbation spectrum and its non-Gaussianity is an interesting 
problem in inflationary cosmology, given our total ignorance about the initial conditions. The 
first steps towards understanding the dynamics of such a scenario with a generic Bianchi I 
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metric were undertaken in [8, 6] and the power spectrum was numerically determined. In 
[10], analytical expressions for the power spectrum of the scalar and tensor perturbations, in 
the case of an axially symmetric Bianchi I metric at early times, were derived for a regime of 
large momentum modes. 

As discussed in [21, 10], a particular axially symmetric Bianchi I background geometry (de- 
scribed in the next section) admits a WKB solution for perturbations at early times for modes 
in a certain high-momentum regime. On matching the WKB solution with the solution in 
the de Sitter phase, one can describe the late time dynamics of the curvature perturbation in 
terms of a non-standard ground state (essentially an excited state on the BD vacuum). Recent 
work [15, 14] claims that particular non-BD vacua can generate visible f^f' 1 for single-field 
inflation. The computation carried out in this work gives an example of a non-BD vacua 
for which this enhancement does not happen. Non-Gaussianity in the scalar perturbation 
spectrum for the non-planar modes in the high-momentum regime (k x ~ k y ~ k z » Hj) were 
studied in [21]. 

In this work, we primarily study non-Gaussianity in the scalar perturbation spectrum for 
another interesting high-momentum regime, viz. the planar modes where k y ~ k z 3> k x , 
where the y and the z axis define the plane of axial symmetry of the particular Bianchi I 
metric under study. We also analyze non-Gaussianities from three-point functions involving 
momenta in different regimes, viz. planar and non-planar. In addition, we study the power 
spectrum for modes in the planar regime and present an explicit analytical formula which 
explains several behaviors found numerically in [6] 

The paper is organized as follows. In the next section, we discuss the anisotropic back- 
ground geometry in which the inflaton evolves and study the classical equations of motion. 
Section 3 deals with the study of cosmological perturbations in the high-momentum regime, 
the WKB solution and computation of the power spectrum in the two regimes mentioned. In 
section 4, we review the results for the bi-spectrum and the /jvl parameter in the non-planar 
regime and in section 5 calculate the analogous quantities for the planar regime. We discuss 
the case of "mixed" modes (i.e. correlation functions involving planar as well as non-planar 
modes) in section 6 and present the conclusion in section 7. 

2. The Model 

We consider a theory of Einstein gravity with a minimally coupled single scalar field given by 
the following action, 



where the background metric is chosen to be an axially symmetric version of the Bianchi I 
metric: 



S 




ds 2 = -dt 2 + e 2p {dxf + e 2 \dy 2 + dz 2 ) 



(2.2) 
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In contrast to the FRW case where one has a single Hubble constant, we have two Hubble 
constants, which we choose to define as follows: 



H=fi±^ t h = ^l (2.3) 

The classical dynamics of the system specified by the action (2.1) constitutes a strongly 
anisotropic expansion at early times (parametrized by h) followed by eventual isotropization 
at a time-scale t « ti so = For t S> U so , the universe enters a phase of de Sitter expansion. 

Note that h, which, roughly speaking, is a measure of the rate of anisotropic expansion 
vanishes in the isotropic limit (p = $) so that we are left with a single Hubble constant. 

In terms of H and h, the independent Einstein's equation and the equation of motion for 
the scalar field reduce to the following set of equations: 

H + 3H 2 = V(4>) (2.4) 
3H 2 - h 2 = ^<j> 2 + V(4>) (2.5) 
(t> + 3Hj> + V'(<j>) = (2.6) 

The time-evolution of h can easily be derived from the above equations, 

h(h + 3Hh) = (2.7) 

In the anisotropic phase, h ^ 0, which leads to the equation of motion, 

h + 3Hh = (2.8) 

For a general V(4>), one can only obtain approximate solutions to the above system of 
equations. However, for a constant V, one can exactly solve the coupled differential equations 
for H, h and <\> as follows, 

H = coth(v / 3~V 7 t) = Hi coth(v / 3yt) 

h = ±W \=- (2.9) 

sinh( V / 3yt) 

= 

In the above solution, the constants have been chosen such that the metric approaches a 
Kasner solution in the limit t — > + . The ± sign in the solution of h indicates two different 
branches in the solution space (distinguished, among other things, by their behavior in the 
Kasner limit). It turns out that only for the positive branch, one can impose initial con- 
ditions on the cosmological perturbations at early times via the usual WKB approximation 
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[11]. Hence, in this note, we will focus exclusively on this class of backgrounds. 



Now, for a given non-trivial V ((/)), the slow-roll condition ((f) « 0) will imply that at early 
times V(<p) is nearly constant with time, provided H(p 2 — > at early times. This condition is 
obeyed by all common inflaton potentials and hence the above solution (2.9) can be trusted 



for a non-constant potential in the t — > + limit. As an example, consider V(4>) 
which H and 4> have the following asymptotic forms at early times, 

m 2 (\>\t 2 



for 



H = 



1 

3t 



1 + 



+ 0(m 4 t 4 ) 



1 - 



m?t 2 



+ 0(mH A ) 



In this case, H<j) 2 ps t — > 0, so that V is essentially constant at early times. 



In this Kasner limit, the metric reduces to the following form, 



ds 2 Kasner = -dt 2 + (Wt) 2 (dx) 2 + {dy 2 + dz 2 ) 



with p = \, /3 = 0. 



(2.10) 
(2.11) 



(2.12) 



This is the gravitational background in which the cosmological perturbations evolve at 
early times. The solutions for the background equations of motion suggest that the universe 
starts life with a very strong anisotropy (h — > \ at early times) which is smoothed out very 
fast by the inflaton potential. The universe then enters a phase of usual isotropic inflation. 



3. Power Spectrum: The Planar and the Non-planar regime 

The computation of the spectrum for cosmological perturbations for a generic anisotropic 
background has two significant differences with the corresponding computation in the isotropic 
case: 

(1) The existence of a WKB solution for the modes of a given wavelength at asymptotically 
early times (t — > + ) is not guaranteed , since for certain backgrounds any mode may become 
super-Hubble [8] in this limit. 

(2) The 50(3) scalar and tensor perturbations are coupled for a generic wavelength at times 
t < t iso [6] . 

As commented in the previous subsection, the particular choice of the "positive branch" 
background solves (1). The positive branch metric in the Kasner limit is, in fact, a patch 
of Minkowski space-time [9] and this can be seen as follows: Let u = t smh(Vx) and v = 
t cosh(yx), so that, 

ds 2 = -dt 2 + (Wt) 2 (dx) 2 + (dy 2 + dz 2 ) = -dv 2 + du 2 + (dy 2 + dz 2 ) (3.1) 
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Therefore, this particular background could admit a WKB solution for perturbations at early 
times. Note that coordinate invariants such as the curvature and the Weyl tensor are time 
independent and non vanishing and hence the space is never Minkowski. Thus, when com- 
puting quantities that depend on derivatives of the metric one should be careful if using (3.1). 

In [6], the perturbations in the anisotropic phase were parametrized in terms of the 
variables V 7 , H x , H + (which reduce to the usual gauge-invariant variables in the isotropic 
limit, V becomes the Mukhanov variable and H x , H + become the two polarizations of the 
tensor modes) with the following equations of motion: 

H';+u 2 x H x =0 (3.2) 

(3.3) 





where the derivative is with respect to the conformal time r\ and the frequencies of the coupled 
system are given as, 

"ii = e 2 »(p 2 x +p 2 y -2pP + ---) (3.4) 
<4i = e 2 "(p 2 x +p 2 y -2p$ + ...) (3.5) 

\2PpI + {p + P)pI) \ M p ) 

The pi are the physical momenta, p x = kie~ 2f) ,p y = k y e~ 2 ^ , where we have chosen the y- 
axis to coincide with the projection of the wavenumber /momentum vector on the plane of 
axial symmetry. The ellipsis in the above equations indicates terms subleading in the limit 
t — > + . To the leading order at asymptotically early times, uj 2 -^ e~ 2p ,uj 2 2 e~ 2p « 1/t 2 while 
u 2 2 e~ 2p « t 2 and as a result, the mixing terms can be neglected. Therefore, in this limit, 
the scalar and tensor perturbations decouple as in the isotropic case. This feature is not 
surprising since, as seen earlier, the metric is that of a flat space-time. 
The frequency of the decoupled H x mode is given as. 



loI = e 2p e^ 



(Pi +P 2 y) 2 



(3.7) 



Now, recall that in the isotropic case, the equation of motion of the Mukhanov variable v 
is identical to that of a scalar field /(after a redefinition / — > a(rj)f) evolving in the same 
background, as long as the slow-roll conditions are obeyed (as a result of which ^ « lt(r)p 
with derivatives in terms of conformal time). In other words, the classical solutions of the 
scalar field / and the curvature perturbation £ are related as / = in the slow-roll regime. 
Therefore, in all quantities of interest, one can substitute curvature perturbation by a solution 
of the scalar field equation of motion, up to a well-defined normalization (specified above). In 
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the anisotropic case, the scalar mode is given as V = exp2B\5d> H = — irr^ — 7^ilA in terms of 

the scalar field and metric fluctuation [6]. Since V approaches the usual Mukhanov variable 
in the isotropic limit, the above argument is also true for the scalar mode evolving in the 
anisotropic background in the limit t 3> ij so when the universe enters a late-time de Sitter 
phase. 

At asymptotically early times, V reduces to purely a fluctuation in the inflaton field in the 
Kasner background. In the limit t — > + , . 2 - p ^ 2 — ^ « t 2 , so that V ~ e 2 ^b<\> (with 

PPy~\~P{^Px~^Py) 

e 2 @ ~ constant). If / denotes a scalar field evolving in the Kasner background, then V sat- 
isfies the same classical equation of motion as e( p+2 ^/ 3 / at early times as long as slow-roll 
conditions hold. 

Thus, both at early and late times, the curvature perturbation can be understood in terms 
of this scalar field /, evolving in the background given by the axially symmetric Bianchi I 
metric, which behaves as a Kasner metric at early time and de Sitter at very late times . 
Therefore, if one can find a WKB solution at early times, one can obtain an approximate 
classical solution at late times by a standard matching procedure at some intermediate time. 
This approximate classical solution will specify a particular vacuum state of the curvature 
perturbation field at late times and can then be used to compute correlation functions in 
the "in-in" formalism. It is important to note that, in this scheme, the entire information of 
anisotropy is encoded in the vacuum state of the theory. 

We are interested in the solution of the above set of coupled equations for the perturba- 
tions at early times in the regime of high momentum : |k| 3> Hj, which essentially implies 
that the modes are deep inside the horizon at early times and exit the horizon at times com- 
parable to the time of isotropization or later (the case of planar modes is slightly more subtle, 
as explained later). 

As discussed in [10], the high-momentum modes can be divided into two distinct cate- 
gories depending on features of their early-time solutions: 

(1) Non-planar modes i.e. k x « k y « k z S> Hj : A WKB approximation is valid at early 
times for this range of high-momentum modes and therefore the late time de Sitter solution 
for the perturbations can be obtained by matching this WKB solution at an intermediate 
time t* = 

(2) Planar modes i.e. k y ,k z 3> k x : The WKB condition is violated at very early times and 
restored at later times. In this case, one first finds a solution to the equation of motion in 
the regime where the WKB condition is violated and matches it to the WKB solution at a 
certain time t x (see Appendix B). This WKB solution is then matched with the late time de 
Sitter solution at t = t* in a way similar to the non-planar case. 

Note that in the non-planar case, the condition |k| S> Hj implies that all the modes are 
sub-horizon at sufficiently early times. On the other hand, for planar modes with k x ~ 0, 
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the modes will be super-horizon in the x-direction no matter how early in time one goes 
back. For this regime of wave-numbers, one cannot impose initial conditions on the modes 
in a consistent way in our formalism. In fact, one can demonstrate that perturbation theory 
breaks down well before you reach this range of wave-numbers. In subsection 3.2, we will 
derive a rough lower bound for k x /Hi imposed by the validity of perturbation theory. 



3.1 Power Spectrum: Non-planar Mode 

In this subsection, we review the computation for the non-planar modes (for details, see [21]). 
As explained in Appendix A, the modes in this regime have a WKB solution at early times, 
which is used to obtain the precise solution in the de Sitter phase by matching at a time 
t ~ 

b * ~ V Hi ■ 

In the de Sitter phase, the classical solution to a scalar field / is given as, 



'2k 3 



'2k 3 



(3.8) 



where 4>±(rj) = (krj ± i) exp (±ikrj), with rj being the usual conformal time as defined in a de 
Sitter universe. 

The coefficients A + and A- are given as, 



A^ 



A. 



(2 - e 2 ) 
2 



• I 1 ) • ()(:'• 



e 3 + 0(e 4 ) 



exp 



exp 



(3.9) 



where e = \J~^~ and r = y^jV+jjJ _ 

The WKB approximation is valid when e <C 1 and we have retained terms up to order e 3 , 
which is the minimal order at which any signature of anisotropy appears. 
^From the above solution, it immediately follows that, 



l/l 



2 



H"2 tj-2 

f±L\A^-A \ 2 -^L 

2k 3\ A + A ~\ ~ 2k 3 



3/2 



cos 2 




(3.10) 



Since the curvature perturbation is related to the scalar field / as / = the two-point 
correlation function of the curvature perturbation is given as, 



<Ck(*)Ck'(*)> ~ (2vr) 3 ( 5 3 (k + k / ) 



2k 3 



1+I3T 




= (2vr) 3 5 3 (k + k' 



^ F(k, cos 6) 
~2k~ 3 



(3.11) 
(3.12) 
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where 



p{t e )e p ^^k (3.13) 
cos# = ^ = \A ~r 2 (3.14) 



~k 

p(t e ) ~ Hi (3.15) 



The power spectrum obtained here is a slightly general form of [20] and admits the parametriza- 
tion, 

P(k) = P(fc) (l + + 9(k)(k.n) 2 ) (3.16) 

where P(k)o is the usual nearly scale-invariant contribution while n is an unit vector in the 
direction which breaks the rotational invariance. 

The spectral index for the curvature perturbation is then given as, 

n a = fc-^-log[F(fc,cos0)] « -^log[F(k,co S 9)} 
dk p e dt e 

« 2(r? - 3e) + Q - cos 2 0^ sin (2e Pe/2 ) + 0(e~ 3pe/2 ) 

(3.17) 

Note that, e _pe ~ e _ * e//tiso . One can easily estimate the magnitude of the correction term 
arising purely due to early-time anisotropy. 

Here, the time t* (time at which we match the WKB solution with the de Sitter solution) 
obeys t* > ti so , such that e p ('*) » 1. In addition, it was shown in [21] (also explained in 
Appendix A) that t e « 2t*. Therefore, for t* = 5^ so , for example, we have t e = 10U SO , which 
implies that the correction term is of the order e~ Pe ~ e -*e/*i SO _ e -io ^ iq~ 5 , while the 
slow-roll terms are of the order 10~ 2 . In this case, we are looking at a regime of momenta 
where £ » 10 5 or, w 10 _1 M P (il e ~ 10~ 6 M P ). 

The contribution of anisotropy to the two-point function could be larger for the regime of 
lower momenta (for which the WKB approximation still holds). For t* = 2ti so , we have 
t e = Ati so which gives a correction term of the order of e~ te ^ tl3 ° = e -4 ~ 1/50 - comparable 
to the slow-roll terms. In this 50, such that the WKB condition is still obeyed. 

However, since these modes have much larger wavelengths, with k ~ 5 x 10~ 5 M p , they may 
well be outside the present regime of observation. In any case, we will more be interested 
in the momentum regime where the contribution of anisotropy to the two-point function is 
negligible and investigate its possible observable signature in the three-point function. 



3.2 Power Spectrum: Planar Mode 

In this subsection, we study the two-point correlation function of the curvature perturbation 
and the corresponding power spectrum for the planar modes. As mentioned earlier and 
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explained in Appendix B, imposing initial condition on these modes involves solving the 
classical equations of motion at early times when the WKB condition is violated. This 
solution is matched to the WKB solution at later times, which, in turn, is matched to the 
late time de Sitter solution. 

As before, in the de Sitter phase, the classical solution to a scalar field / is given as, 

„ A + (k)#j , . . A_(k)fTj , . , 

/ = -^y^Mv) + u= Mv) (3.18) 

where 4>±(rj) = (krj ± i) exp (±ikr]), with r\ being the usual conformal time as defined in a de 
Sitter universe. 

The coefficients A + and A-, obtained by matching procedure (explained in Appendix B), are 
given as, 

e i(ir/4-bk/H z ) 

A + = i 

■J \ _ e -2it(ask/H I ) 

-7r(osfc/Hj) -i(it/4-bk/Hr) 

A- = = (3.19) 

J\ _ e -27r(asfe/if/) 

where s = k x /k < 1 is a measure of "planarity" of the mode while a = and b = 

22/3 /j{T(l/3) 

— 3 r(5/6) are numerical constants. Note that A± obey the usual normalization condition, 
viz., 

\A+\ 2 - \A_\ 2 = 1 (3.20) 

Given that the classical solutions of / and ( are related as, / = -^Q , the two-point correlation 
function of the curvature perturbation is given as, 



(Ck(iXk'(i))«(2^ 3 (k + kO^|| 



e L 



r , TT , sm(2bk/H I ) 1 

coth (irask Hi) - —j '—jf- 3.21 

smh (irask/Hi) 



where 



p(t e ) e p( * e) ps k (3.22) 
p{t e ) ~ Hi (3.23) 



The power spectrum for the curvature perturbation then reads, 

PW = FW„ (coth (nast/H,) - (3^24) 

where P(k)o = -± is the power spectrum for standard inflation. 

It is useful to compare the above expression to the power spectrum obtained for the non- 
planar case. For a generic value of s, the power spectrum is not a perturbative modification 
of the standard power spectrum and differs non-trivially from the ACW formula (equation 
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3.16). In other words, the effect of anisotropy can be large even in the limit of high momenta 
in the planar regime and can strongly violate the scale-invariance of the spectrum, as we 
discuss below. 

There are two interesting limits in which the power spectrum can be analyzed. Firstly, 
consider the limit s — > (such that nask/Hj <C 1 although k/Hj » 1). In this case, we have, 

It is instructive to compare the above equation to the behavior of the power spectrum found 
numerically in [6] (equation (35)). For small enough s, the power spectrum precisely diverges 
as 1/s, as predicted in [6]. Moreover, for a given k, the range of s for which this behavior is 
observed, shrinks to smaller size as k grows, simply because s has to be progressively smaller 
in order to satisfy the condition nask/Hf <C 1 as k grows. This is also in agreement with the 
results obtained in [6]. 

The oscillatory behavior of the power spectrum (fig 3 of [6]) for small s (i.e. s <C Hj/k) 
can also be analytically seen from the above formula. Smaller the value of s larger are the 
oscillations of the power spectrum. In the limit of large k, i.e. s S> Hj/k, the oscillations in 
the power spectrum die out and the spectral index approaches the scale invariant formula. 

There is a strong bound on the planarity parameter s due to the fact that perturbation 
theory has to be valid at early times and we derive this now. Since curvature perturbation can 
be understood in terms of the scalar field / evolving in the Bianchi I background, perturbation 
theory will be valid if at all times the total energy of the scalar field is much smaller compared 
to the energy of the background, which (given that the metric only depends on time) implies, 

d 3 xT™ » J d 3 xT™ (3.26) 

where X^jj is the background energy density and TjP^ is the energy density due to the scalar 
field /. Note that 



J d 3 xT™ = J d 3 x ( 



Q) 2 f (k, T )f (-k, r) + gVkikjfiK r)/(-k, r) 



where r is the "conformal time" defined in Appendix B and prime denotes derivative with 
respect to r. Now, plugging in the classical solution of / at early times, i.e. /(k, r) = 
exp [— 3iak±T + tip] (also derived in Appendix B), we have, 



d 3 k 



—&iak x T+iip 



(3.27) 
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where ip is a time-independent phase factor: e 1 ^ = e - iak x/H I \n(ak/H I ) _ 

The only way to have a vanishing (or small) contribution from the above integral is to have 
a large exponent for the exponential factor so that integration over k x involves major cancel- 
lations. We therefore require, \k x r\ 3> 1. Since r « ^ In Hjt, this implies, | jyy lni//i| 3> 1. 
Taking i7j = 10~ 6 M P and the minimum time-scale at which equations of General Relativity 
still holds as t m i n = 10M" 1 , we obtain jj- S> 1/10. Therefore, the value of s for which our 
formalism is valid must obey the following (/c-dependent) lower bound: 

j| > 0(1) (3.28) 

The range of wave-numbers for which the power spectrum gives al/s behavior must satisfy 
^ < 1 and is therefore outside the physically sensible regime of wave-numbers. 
The second interesting limit is tt ask /Hi » 1 (with s < 1), such that, 



P(k) » P(k) (3.29) 
where the isotropic behavior of the power spectrum is retrieved. 

Finally, the spectral index for the curvature perturbation can be parametrized as n s = 2(rj — 
3e) + An s , where An s can be derived from the above equation as follows: 

k .7rassinh {irask/Hj) — 26 cos (2bk/Hj) , . , .„ , n 

An s = — -i V / ~, \ ~ kcls coth (irask HA) (3.30 

Hi cosh (yrasfc/^/) - sin {2bk/H I ) v 1 1JJ v J 

For s > Hi Ik such that sinh ask /Hi) ~ cosh (it ask /Hi) ~ e nask / H i ) we have 

An s w -JL e ~™sk/Hi cos (26^;/^) (3.31) 

Equation (3.31) directly gives the regime of planar modes which preserve the scale invariance 
of the curvature perturbation spectrum, 

^ask/Hj > k 

Hi 

For a generic value of s in this regime, the analytical formula for the spectral index suggests 
that the effect of the initial anisotropy might be substantial. However, in order to determine 
whether this effect will be observable, one needs to study the resulting statistical correlation 
between different multipoles of the CMB temperature anisotropics. In the next subsection, 
we numerically calculate the signal-to-noise ratio for individual multipoles for a physically 
reasonable range of values of Hi as an estimate of how large this effect can be. This analysis 
shows that the signal-to- noise ratio is pretty low even for the maximum allowed value for Hj. 
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3.3 Angular Correlation Matrix and Signal- To-Noise Ratio: Planar Modes 

Our starting point is the expression for the power spectrum obtained in the last subsection 

(Cfc(*)Cfc(0) = (2vr) 3 ( 5 3 (k + k / )P(fc)o (coth(7rafc|cos0|) - s " 1 ( 2 ^/^) 
N ' \ smh(7ra/c | cos 6*1) 

= (27r) 3 5 3 (k + k')P(M) (3.33) 

where P(k)o is the isotropic LCDM power spectrum and is the angle between k and the 
anisotropic direction of the Bianchi I metric. 

We now want to connect this with the observable temperature anisotropics in the cosmic 
microwave background radiation. The full covariance matrix for the temperature anisotropies 
is given by 

// ST ST \ 

dQdtt' /— (n,r ? o,xo)^ r (n',r ? o,x )\y^(n)Yg, m ,(n / ), (3.34) 

where the temperature perturbation ST/T(h, rjo, xq) = m a£ m Yi m (h) in the direction n 
observed from xq is related to the power spectrum of metric perturbations through 



ST f r/ 3 k 

— (n, % ,xo) = J ^3C(k,r / ,)A(A ; ,k.n,r ?0 ) e lk - x o (3 . 35) 



The transfer function A(fe, k-n, 770) in fact depends only on the magnitude of the wavenumber 
k = |k| and on the angle between k and the line of sight n, since the initial conditions for the 
evolution of perturbations after the onset of inflation can be considered isotropic [34]. Using 
the line-of-sight approach [26] and omitting 770 from now on for simplicity, we can expand in 
spherical harmonics 

A(fc, k • n) = J^H) L (2L + l)P L (k • n)A L (fc) 

L 

A L 

PL ^' h) = T2LTT) £ Yl M (M)YLM{n) (3.36) 

^ ' M=-L 

and derive a generic expression for an anisotropic initial curvature power spectrum 
— 3P(k)(4vr) 2 (-i) L - L '&L(k)AUk)YZ M (k)Y L , M ,(k) 

^ n ' LL'MM' 



x J dndn'Y; m (h)Y e , m ,(h f )Y LM (h)Yi IM/ (h^ 

= i-^'l J d 3 kP(k)A e (k)A* ei (k)Y; m (k)Y Pml (^)- (3.37) 

To simplify this expression for the axially symmetric power spectrum in our model, we use 
polar coordinates in momentum space and rotate the system so that 6 is the angle between k 
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and z (which now points in the anisotropic direction of the metric) . Performing the (j) integral 
and using Y tm (e,<f>) = ^ J^$^ PP (cos 9) e im <t> we get 

1 

Cu'mm' = 2 6 mm > (-if~ l 'f k 2 dkA e (k)A* e , (k) J d(cos 9)Y lm (9, = 0) Y Vm (9, = 0) P(k, 0)(3.38) 

-l 

To calculate the signal-to-noise, we separate the isotropic and anisotropic contributions 

Chnl'm' = ( a e m a e'm>) = ^mm' Cj> tS °^ + <5Q m £' m ' , (3.39) 

where cf s °^ is the isotropic covariance matrix recovered by substituting P(k, 9) = P (k) 

cf o) = 5 u ,5 mm ,C^ m , = k 2 dkP (k) \A £ (k)\ 2 (3.40) 

and SCimi'm' is the deviation as a result of the anisotropic Bianchi I metric. We start from 
the likelihood function for an isotropic LCDM covariance matrix [27, 28, 29, 30] for a system 
of n Gaussian fields satisfying (xiXj) = c\™°^ 

6 iso ^ = 1 exp f- 1 x T (C , ( is °))~ 1 xY (3.41) 

Given our perturbation C = + 5C we follow [30, 31, 32] to define 

X 2 (C) = -2 [log(£(C)) - log(£( is °))] = x T (c' 1 - (C^)" 1 ) x + log (det C/ det C (iso) ) (.3.42) 
The signal-to-noise for detecting the induced anisotropies is given by the mean value of \ 2 

£)' = <C)> =/d^(C)£<M = T,- ( W - 1) + .og (detC/de,C««>) ,3.43, 

Expanding this to leading order in SC, we can examine the behavior for different multi- 
poles £ 

!)'-?(£)!• (3 - 44 » 



where (given the axial symmetry in our model) 

i'=2 x 7 m=-l' °c u <" 



In Fig. 1 we plot the signal-to- noise as a function of multipole t. 12 The power spectrum 

correction affects small wave numbers k, yielding a signal in the largest scales, or smallest 

multipoles £, for which the statistics is more limited by cosmic variance. We see that even 

for the maximum value allowed for Hj, the signal-to- noise is quite low and conclude that this 

feature in the power spectrum would most likely remain beyond experimental reach. 

To avoid the breakdown of the perturbative treatment in the limit -nak cos 6 /Hi <§C 1, we cut off the power 
spectrum correction in (3.33) at cos# = 5Hi/(ivak) (verifying numerically that the results do not depend on 
this specific choice). 

2 We used transfer functions from CMBFAST [26] using CMBEASY [33] in the calculation. 
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Figure 1: Left: The ideal signal-to-noise per multipole I for different values of Hj between 10~ 5 M p i 
(yielding almost no signal) and 10 _4 M p i (peaking at I <~ 7 as a result of the tradeoff between cosmic 
variance and the signal dominance in large scales). Right: The total signal-to- noise in £ < 30 for 
different values of Hj. The upper bound Hj ~ 10 _4 M pl yields a signal-to-noise ratio of <~ 3.5 and will 
be challenging to detect in a real experiment. 



4. Bispectrum: Non-planar modes 

In this section, we briefly review the results for the 3-point function of the curvature pertur- 
bation and the f^L parameter for non-planar modes (for details, see [21] ). 
For computing the 3-point function for the curvature perturbation, following [2], we consider 
a local (in time) non- linear field redefinition of £: 



2# Sc ' 8^ 2Sc ' V 



C = Cc + "He 2 + h>£ + ^d~H( c d\ c ) (4.1) 



Evidently, this redefinition does not change the quadratic action which implies that Q c and 
£ have the same equation of motion and hence the same classical solution given by equation 
(3.1). Also, since local redefinitions do not yield any enhancement of the 3-point function, it 
is sufficient to compute the correlation function in terms of the redefined field Q c . In terms of 
Cc the leading term (in slow-roll parameter) in the interaction Hamiltonian will be given as, 

Hi = -J d 3 xd V e 3 e (j^J f>C' c 2 d~ 2 C' c (4.2) 

where the prime denotes derivative w.r.t. the conformal time rj (defined in the de Sitter 
phase), and the partial indicated space derivatives. 

We can now use the "in-in" formalism to compute the tree-level contributions to the 3-point 
function, corresponding to the de Sitter vacuum specified by equation 3.9. 
using the usual Feynman rules in the momentum space, the three-point correlation function 
at a conformal time r\ is given as, 

<Cc(ki, r/)C c (k 2 , r/)C c (k 3 , v)) R/L ~ <5 (3) ( k i + k 2 + k 3 )A i? / L (k 1 , k 2 , k 3 , r?) (4.3) 
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where A R and A L are given as, 

A*(ki,k 2 ,k3,»7) = i P drft 




i=i 



A^(ki,k 2 ,k8,»/) = (^(ki,k 2 ,k 3 ,77)r 



(4.4) 
(4.5) 



where Gk ; (?7, 77') = Cd(ki, r7)Q(ki, rf), with CcK^i, f?) being the classical solution for curvature 
perturbation. A R (ki, k 2 , k3, rf) denotes the contribution to the 3-point function correspond- 
ing to the tree- level Feynman diagram with a "right" vertex while A R {k\, k 2 , k3, rf) denotes 
the contribution corresponding to the tree-level Feynman diagram with a "left" vertex. The 
final result for the 3-point correlation function of the curvature perturbation is given as, 

(Cc(ki^)Cc(k 2 , v)Uks, V)) ~ <5 (3) (ki + k 2 + k 3 )[A R (k 1 , k 2 , k 3 , rj) + ^(ki, k 2 , k 3 , r,)] (4.6) 

In the definition of A R (\&i, k 2 , k 3 , 77), one needs to make a choice of r?o - which in standard 
inflationary scenario is taken to be —00. However, we will take 770 to be of the order of the 
isotropization time-scale when the universe has entered an essentially de Sitter phase, thus 
rjo is near the onset of inflation and parallels the choice made in [5] . In the computation 
of late-time correlation functions, the parameters p and (j> can therefore be assigned their 
respective de-Sitter values, which remain nearly constant during inflation. 
In the limit of late times (77 — >• 0), we have, 



A R (k 1 ,k 2 ,k 3 ) 



2 



&=±li=l 



(4.7) 



where the sum extends over all 8 possible linear combinations and F^ i= _\{ki) = ^(|^4+| 2 — 

AlA\*) and F^ =1 {kf) = g(|^_| 2 - A\Al*). 
Therefore, on completing the rf integration, we have, 



A*(ki,k 2 ,k3) = - 



k 2 h 2 



E (II 



&=±i i=i 



where the functions F±(kf) are given as, 



6 



Ft 



4 

2p 



= tp 
1 + 



1 



+ cos 2 9i I ef exp 



1 



1 



+ cos 2 6i ) ef exp 



-2i 

£i 

2i 

Si 



(4.8) 

(4.9) 
(4.10) 



There are two interesting configurations for which this amplitude is enhanced, 

(1) Flattened triangle limit fc 2 ~ fc 3 ~ k±/2 ~ k, such that k\ = k2 + A; 3 . As discussed in 

[21], the contribution of early time anisotropy to the f^L parameter is given as, 



\Af NL \ =e 



(4.11) 
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which, in addition to the standard slow-roll factor, is suppressed by powers of ^, leading to 
an extremely small change in the value for Jnl- 

(2) Squeezed triangle limit k$ <C k\ ps hi ~ fc but |fci — fol 7^ £3. In this case, the change 
in the /nl parameter in the standard case is given as [21], 

W«l -.(*)'£ (W 

which shows that the enhancement for the "squeezed triangle" limit is greater compared to 
the "flattened triangle" by a factor of k/k%. However, since {k/kz) ma x ~ 200 given the range 
of scales measurable in the CMB, |A/^vl| is too small to be observed even for k/Hj ~ O(50). 
For the range of wavelengths for which the contribution of anisotropy to the spectral index 
is negligible, the contribution is obviously far smaller. 



5. Bispectrum: Planar modes 

In this section, we derive expressions for the Jnl parameters for planar modes and study 
these in the regime s > Hj/k , which includes the extreme limit sk/Hi 3> 1, where the 
scale-invariance of the power spectrum is retrieved. 

The 3-point correlation function for the planar modes can be computed similarly in the "in- 
in" formalism using the interaction Hamiltonian given in equation (4.2) corresponding to 
the Bogoliubov transformed ground state specified by equation (3.19). Following the same 
procedure illustrated for the non-planar modes in the previous section, the "right" amplitude 
can be written as, 

^( kl ,k 2 ,k 3 ) = -i^M [°d v > £ fl^WF^kt) (5.1) 



&=±li=l 



where the sum extends over all 8 possible linear combinations ^2iiki and F^ i= _i(fej) 
Sd^VI 2 - AlA\*) and F^ =1 (k) = §(|A!_| 2 - A\Al_*). 



Therefore, on completing the rj integration, we have, 



^( kl ,k 2 ,k 3 )=-^s y: (u^^TTki 1 -^ ^) 

*r £i=±l i=l ^ l 1 



(5.2) 



2fc 3 

Let us define F^.(ki) = -j^-F^(ki), which can be written as explicit functions of wavenumber 
using equation (3.19) as follows: 

M _ e -Tv(ask i /H I ) e -i(iv/2-2bki/H I )-\ ^aski/Hj _|_ ^gi26fei/if/ 
F ~^ = i _ e -2n(ask t /Hj) 2 sinh (vrasfc/tf/) ^ 

F+(ki) = - e -™'°i/ H i e iW2-™°i/ H i)F-(k i ) (5.4) 
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As for the non-planar case, the amplitude computed above is enhanced for the flattened tri- 
angle configuration as well as the squeezed triangle configuration and are analyzed as follows: 



Flattened Triangle Configuration 

In this case, the enhancement appears when the wave-numbers satisfy ]P = 0) so that 
the exponent of the exponential term in equation (5.2) vanishes. 

We choose £12 ~ &3 ~ fci/2, setting k\ = &2 + k%. This implies that the set {&} contributing 
to the enhanced bispectrum can have values (1,-1,-1) and (—1, 1, 1). Therefore 

A*(k,k,k) = -^tF3(k) Vo (e-™ sk / H 'e- i2bk / H ' - ie~ 2 ™ sk / H ' e - mk ' H ') (5.5) 
8k 5 (f) 2 . 



Using the definition of /j^x*: 



fg? = (A R (k,k,k) + c.c.)/p(ky 

we have 



A .flat = 4>l ( cos (2bk/H T ) cosh (Trasfc/ffj) - sin (4&fc/ffj) i 
!nl ' Hf V (cosh {-Kask/Hi) -sin {2bk/H I )f 1 



In the s > Hj/k limit, 



A ff lat ~ ^ -(nask/Hj) 
J NL ~ jj2 e 



\MUl\ « 



L i 

which implies that the f^L parameter is suppressed by the exponential factor in addition 
to the slow-roll parameter. In the regime where the power spectrum is scale invariant, 

e (nask/Hj) > k j H ^ we have 

f fiat l ^ 4>l H i 
Hf k 

This configuration, therefore, does not produce any appreciable enhancement in the /jvl pa- 
rameter in any regime. 

Squeezed Triangle Configuration 

For the squeezed triangle configuration, ks <C k\ « ki ~ k and the corresponding amplitude 
is given as, 

where the sum is over the set £j : (1, —1, —1), (—1, 1, 1), (1, —1, 1), (—1, 1, —1). 

Let S3 denote the planarity parameter for the smaller vector k3, while s denotes the planarity 

parameter of the vector k. Using the definition of f^if- 

sqzd = A R (k,k,k 3 ) 
Inl P(k)P(k 3 ) 
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we have 

A./ 



sqzd = f <K_ \ ( k_ \ 1 

NL \Hf ) \k 3 J smhirask/Hjicoshirask/H! - sin 26fc/fl» 



a(fcs, r/o, s 3 ) + e™ sfc /^ fc 3 , a 3 ) + e-* aak ' Hl -y(k, h, r? , a 3 ) 

(5.8) 



(cosh nas^kz/Hj — sin2&£: 3 /.f/j) . 
where the functions a(k 3 , rjo, s) , /3(k, fc 3 , rjo, s 3 ) and ^(k, k 3 , rjo, S3) are given as 



a(h,Vo,s 3 ) = -2( e 7raS3k3 +ie l2bk3/Hl ) \{1 - jwk^ + ie -™ S3 k 3 /Hi e -i2bk 3 /H T ^ _ e -i W fc 3 ) 
/3(fc, fc 3 , %, s 3 ) = z(e™ S3fc3 +ie i2Wc3 /^) [(1 - e^ 3 ) + ie -^s 3 k 3 /Hj e -i2bk 3 /Hj ^ _ e -i Vo k 3 ^ 



e i2bk/H I+c _ a 



+ C.C. 

(5.9) 
(5.10) 

l(k, k 3 , Vo, s 3 ) = -i(e waS3k3 +ie i2bk3/Hl ) \(1 - e ir >° k3 ) + ie -^sk 3 /Hj e -i2bk 3 /Hj ^ _ e ~i V0 k 3 

(5.11) 

As before, in the regime s > Hj/k, one can estimate the order of magnitude of the f^L 
parameter: 

\MTl\ « (§2) (|) *-™ Sk/Hl (5-12) 



Evidently, the squeezed triangle configuration can lead to larger non-Gaussianity in this 

< k ' 

^ fc 3 ■ 



regime compared to the flattened triangle case. However, since (ir~ )|max ^ 

0(100), the 



enhanced /nl can at best be of O(10 



6. Bispectrum: Mixed case 



-2\ 



We now consider the case where some of the momenta {fej} belong to the non-planar regime 
while others are in the planar regime. As before, we will only compute the leading contribution 
to /nl due to the background anisotropy. 

The amplitude from the "right" vertex in this case can be written as, 

^( kl ,k 2 ,k 3 ) = -?^p £ (f[GM)r^(l-e*»Xi**) (6.1) 

where G^(ki) = F^(k{), if fcj is in the planar regime, and G^(ki) = F^ p (ki), if ki is in the 
non-planar regime. Recall that the functions F^.(ki) are given as 

F^ih) = ^- 3 (\A\\ 2 - AIA\*) (6.2) 
F^ =1 (h) = ^- 3 (\Al\ 2 - A\Al*) (6.3) 
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For the planar case, the functions A± are obtained from equation (3.19) , while for the non- 
planar case these are given by (3.9). Also, the functions F^(ki) are defined as F^.(ki) = 



H - 

The bispectrum has enhancements for both flattened and squeezed triangle configurations, 
which we catalogue below. 

Flattened Triangle Configurations 

(1) One non-planar momentum 

We choose ki = k2 + k3 for the flattened triangle configuration. There are two different cases 
to consider: 

(a) ki and k2 are planar and k% is non-planar 

In this case £i = 1,^2 = —1,^3 = —1- The "right" amplitude then reduces to 

A*(k,k,k) = ^Fl(k)F p (k)F^(k)^k V0 ) (6.4) 
Therefore, the change in the Jml parameter is given as, 

flat _ A R (k,k,k) + C.C. 

P^(k) + 2P p (k)P np (k) 

(0g \ sinh irask / Hi cos 2bk / Hi 

H^ J (coshnask/ Hi — sin26/c / Hi)(2 smhirask / Hi + coshirask/ Hi — sin2bk/Hi) 

(6.5) 

(b) k2 and k3 are planar and ki is non-planar 

In this case £i = — 1, £2 = 1, £3 = 1- The "right" amplitude then reduces to 

A*(k,k,k) = ^-F" p (k)^(k)^(k)(^ ) (6.6) 



Therefore, the change in the f^L parameter is given as, 

a .flat A R (k,k,k) + C.C. 

Jnl Pi(k) + 2P p (k)P np (k) 

f <j> 2 e \ - sin \bk /Hi + 2e~™ sk / Hl cos 2bk /Hi 



I H"j J (coshnask/Hi — sin2bk / Hi)(2 sinhirask / Hi + coslnrask/ Hi — s'm2bk/Hi) 

(6.7) 

The total change in f^L is given as the sum of the contributions (a) and (b). In the limit 
s > Hi /k we obtain 

12 



A ff lat ~ ^ e ~^asklHi 
^JNL ~ T72 e 



Hi 
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so that no appreciable enhancement can be observed for this configuration. 



(2) Two non-planar momenta 

We again choose ki = k2 + k3 for the flattened triangle configuration. In this case, ki is 
planar while k2,k3 are non-planar corresponding to £1 = 1,^2 = — 1, ^3 = — 1. 
The "right" amplitude in this case is given by 

A*(k,k,k) = ^^(k)^(k)^(k)(ife») (6.8) 

Therefore, the change in the f^L parameter is given as, 



flat _ A R {k,k,k) + C.C. 

P|(k) + 2P p (k)P np (k) 



A2 



cos2bk/H I 



Hj J (sinhirask/Hj + 2coshirask/Hj — 2sin2bk/Hj) 



(6.9) 



For the same reasons as explained in case (1), non-Gaussianity from this configuration is also 
severely suppressed. 

Squeezed Triangle Configurations 

(1) One non-planar momentum 

Taking ki = k2 + k3, with k% <C k±, k<i we again have two cases to consider: 

(a) k3 is non-planar : This corresponds to the case where £1 = ±,£2 = —1,^3 = —1, with 
ki,k2 both planar. The "right" amplitude is given as, 

A*(k,k3) « jMt^ 1 ~ e- i%fes )^(k)F p (k)F^(k 3 ) (6.10) 

Therefore, the change in the f^L parameter is given as, 
sqzd _ ^ R (k,k 3 ) + c.c. 



A/ 

~ h (sinhvrasA;/F / )(cosh7rasA;/ J ff / - sin26fc/#j) ^ ~ COS7] ° k3 " ) 

+ e vash l Hl {- sm2bk/H! + sin {2bk/Hr + m k 3 )) - e ~™ sk / H ' ( S m2bk/ Hj + sin (26fc/fTj - 770&3))] 

(6.11) 

(b) k3 is planar : There can be two contributions in this case depending on whether £3 
is +1 or —1: 
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(i) £3 = -1 
In this case, 



A R (k,k 3 ) « " e-*»*»)^(k)^(k)#f (ks) (6.12) 



The change in the /jvl parameter is given as, 



A , sgz d = A R (k,k 3 ) + c.c. 

P p (k)P p (k 3 ) + P np (k)P p (k 3 

ii| J k 3 (coshiras 3 k 3 /Hi — sm2bk 3 /Hi)(smh.Trask/ Hi + coshnask/Hi — sin2bk/Hj) 

(6.13) 



In this case, 



ff 6 

.4 R (k,k 3 ) » ^-^-(1 - e-^ fc «)F" p (k)^(k)^(k 3 ) (6.14) 



The change in the /nl parameter is given as, 



sqzd _ A R (k, k 3 ) + C.C. 

P p (k)P p (k 3 ) + P np (k)P p (k 3 ) 

( 4>l \ k (1 - e _i '» fc3 )(e-' raS3 * 3 / frj - j e -<2Wfc/ff/)( e -™*k/fTj _ ie -i2bk/Hjj + c c 
\Hj I ks (cosh iras 3 k 3 / Hi — s'm2bk 3 /Hi)(s'mhirask/Hi + coshirask/Hi — sin2bk/Hi) 

(6.15) 

The net change in the /at l parameter in the case of one non-planar momentum with k 3 being 
planar is the sum of the contributions (i) and (ii). 

(2) Two non-planar momenta There is only one case of interest in this category, which cor- 
responds to £1 = 1,^2 = — 1,^3 = — 1 - k2,k 3 are non-planar. 

A*(k,k3) « - e-^)Fl(k)F^(k)F!l p (k 3 ) (6.16) 

The change in the f^L parameter is given as, 

f sqzd _ A R (k,k 3 ) + C.C. 



A/. 



x L ~ P np (k)P np (k 3 ) + P p (k)P np (k 3 ) 

( 4>l\ k e -w«*/#i(i_ cog _ gin 2bk/Hi + sin (2bk/Hi + r) k 3 ) 
\Hj J k 3 (sinhirask/Hi + cosh irask/Hi — sin 2bk/Hi) 

(6.17) 
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In the regime sk/Hj > 1, one therefore has the following leading order contribution to the 
Inl parameter for all the above "mixed" cases: 

MTl - (|>) y/™ s k/Hl (6 - 18) 

which is precisely the same as the contribution of the squeezed triangle configuration of 
purely planar modes (5.12). Choosing |% ~ 10~ 2 , ~^\max ~ 200 and sk/Hj ~ 2, we have 

Af s ^ d ~ O(10- 2 ). These two cate gories of squeezed triangle configurations maximize the 
enhancement in the bispectrum due to early time anisotropy. 

7. Conclusion 

This note presents an analytical approach for studying the signature of a strongly anisotropic 
pre-inflationary phase in the spectrum of late time cosmological perturbations for a range of 
high momentum modes which typically exit the horizon after isotropization. Starting with 
the previously known case of non-planar modes (k x ~ k y ~ k z ) [21], we have worked out in 
complete detail the case of the planar modes (k x <C k y ,k z ). While non-planar modes always 
admit a WKB description at early times, WKB condition for planar modes is violated at 
very early times and restored later in course of their evolution. In either case, the late time 
dynamics of the modes is characterized by an excited state (built on the standard Bunch- 
Davis vacuum by a Bogoliubov transformation) which carries all the information about the 
pre-inflationary anisotropic phase. The problem of space-time anisotropy at early times can, 
therefore, be recast into a problem of quasi de Sitter evolution with a non-standard vacuum 
state at late times. We therefore study the correlation functions (two and three point func- 
tions) of curvature perturbation in this excited state and calculate the power spectrum and the 
leading order contributions to the f^i parameter for both planar and non-planar modes due 
to anisotropy. In addition, we compute the leading order contributions to the f^L parame- 
ter in cases where the 3-point correlation function involves both planar and non-planar modes. 

In the power spectrum for non-planar modes, the contribution of anisotropy is essentially 
perturbative in nature, since the corresponding (leading order) correction term in the spectral 
index for these modes is suppressed by a power of e~ te ^ ti3 °. As a result, for t e /ti so ~ O(10), 
for example, the correction due to anisotropy is severely suppressed. However, in the case 
of planar modes, the power spectrum for a generic wavenumber k and parameter s (defined 
as s = k x /k), strongly breaks scale invariance. In the limit s <C Hj/k, the power spectrum 
appears to diverge as P(k) ~ f(k)/s, which agrees with predictions from numerical computa- 
tions in [6]. Also the range of s (for a given k) for which the aforementioned s- dependence of 
the power spectrum is valid shrinks with growing k. However, from considerations of validity 
of perturbation theory, we show that the inequality sk/Hj > 1 must be satisfied and therefore 
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the regime s <C Hj/k is not physically relevant. 

For generic modes in the allowed regime sk/Hj > 1, the analytical formula for the spec- 
tral index suggests that the effect of the initial anisotropy might be substantial. We compute 
the resulting statistical correlation between different multipoles of the CMB temperature 
anisotropies and numerically calculate the signal-to-noise ratio for individual multipoles for 
a physically reasonable range of values of Hj as an estimate of how large this effect can be. 
This analysis shows that the signal-to-noise ratio is pretty low even for the maximum allowed 
value for Hj. 

The bispectrum for non-planar modes computed in the non-standard vacuum state is 
enhanced in the flattened triangle configuration (&2 ~ ks ps k±/2) as well as the squeezed 
triangle configuration (£3 <C k\ « fo), with the latter yielding a larger /jvl parameter. How- 
ever, in both cases, the /tvl parameter is suppressed by factors of e - 3 W 2 <i SO m addition to the 
slow-roll parameter, which constrain the non-Gaussianity in the scalar perturbation spectrum 
to be extremely small. For planar modes, the bispectrum is similarly enhanced for the flat- 
tened triangle as well as the squeezed triangle configuration. In each case we derive the /jvl 
parameter for a generic planar mode of wavenumber k and parameter s and study its behav- 
ior in the regime sk/ Hj > 1 . The squeezed configuration again gives a larger enhancement 
in the spectrum compared to the flat configuration. However, Af^^f « {jfi){-^) e ~ %ask ^ Hl 
involves strong exponential suppression as a result of which it can at best be of O(10 -2 ). 
In the regime of scale-invariant power spectrum i.e. sk/Hj » 1, Af^ d is vanishingly small 
as before, because of larger exponential suppression. The condition of scale invariance, there- 
fore, naturally leads to unobservable non-Gaussianity for the planar modes. 

Finally, we also analyze non-Gaussianity originating from three-point correlation func- 
tions involving both planar and non-planar modes. As before, one expects enhancements from 
the flattened as well as the squeezed triangle configuration with the squeezed triangle config- 
urations producing larger non-Gaussianity. The leading order contribution of the squeezed 
triangle configurations in the "mixed" case to the Jnl parameter turns out to be the same as 
the purely planar case: Af^^f ~ (jjk)j^e~ nask / Hl . For the choice of parameters described 

m the previous section /\ f j^j^ I max 

o(ht 2 ). 

To summarize, our computation suggest that primordial anisotropy will not produce an 
observable f^L in the near future. 
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8. Appendix A: Non-planar Modes 



In this section, we study the evolution of the non-planar modes for a single, massless scalar 
field (referred to as / in the main body of the paper) minimally coupled to gravity in an axially 
symmetric anisotropic space-time, with a positive cosmological constant. In particular, we 
construct a WKB solution for these modes at early times and show how one can match it 
with the general solution at late times to obtain equations (3.9). The treatment essentially 
follows [9]. 

We consider the following action for the scalar field, 

S = - J d^x^-g Q<7^<R0 + V^j , {M 2 p = 1) (8.1) 

where the background metric is chosen to be an axially symmetric version of the Bianchi I 
metric: 

ds 2 = -dt 2 + exp {2p){dxf + exp {2f3)(dy 2 + dz 2 ) (8.2) 
where p, (3 are known functions of time: 

1 3Hjt , TT 

p = -mtanh — - — smh3ii/t 

1 sinh3i?/t 

^ = 3 ln tanhW ^ 

with Hi = ^/f . 

Define p = a — 2/3+ , /3 = a + j3 + and a new "time" coordinate r, analogous to the conformal 
time in the isotropic limit, as, 

dr=% (8.4) 

From equation (8.3) and the definition of the coordinate r, one can derive, e 3a = e p+2/3 = 
sinh3gjt= ginh } 3Hir . 

It can be easily seen that as t varies from + to oo, r varies from — oo to 0~. In this time 
coordinate, the equation of motion for a mode 4>k is given by, 

(j^ + ^) 2 )<^ = (8-5) 

The frequency squared is given as, 

u{r) 2 = 2 ^{l-r 2 x) (8.6) 

where r 2 = ky ^ z and x{t) = 1 — e 6HlT = exp (— 6a)(y / exp (6a) + 1 — 1). Evidently, x{t) 
varies from 1 to as r changes from — oo to 0. 
Equation (8.5) has a WKB solution: 



mKB = —j= exp 

■'2oj j r 



\-i f dr'u] (8.7) 

J TCI 
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where Cj has to be determined from the equation, 



The WKB approximation holds as long as the WKB parameter 



^3 



Hi l-x(r) 



+ 1 < 1 (8.9) 



k (x/2) 1 /3(i_ r 2 x ( T ))i/2 ^i_ r 2 x ( r ) 

The choice of the WKB solution above is obviously equivalent to imposing a particular initial 
condition on the modes of the scalar field at early times. This can be seen directly by 
analyzing the early time behavior of the classical solution for the scalar field. Firstly note 
that Cj ~ u — > k x in the limit r — > — oo (or t —> 0+) and in this limit r and t are related as 
r = jfij l n Therefore, in the early time limit, the time-dependence of the WKB solution 
is given as follows, 



1 r . k x , 3H I t 1 



4>wkb ~ -j== exp[-i— - log — — ] ~ t 3H i (8.10) 
V lk x irij I 

Now consider the equation of motion of the scalar field at early times, 

k + \k + (k 2 y + k 2 z + ^)0 fc = (8.11) 



Define z = ln ^J~k^~+~k^t, so that the equation reduces to 

+ (e 22 + = (8.12) 

In the limit t — >■ 0+, z — )■ — oo, so that the exponential term drops out of the above equation 
and we obtain a solution of the form, 

<f> k (z) = A{k)e~ i ^ z + B{k)e i ^ z (8.13) 

Choosing B(k) = 0, we find that (f)k(z) ~ t 3H i , confirming that the WKB solution has 
the same time-dependence at early times as expected from the classical solution subject to a 
certain initial condition. 

In fact, equation (8.11) has the general solution, 

Mt) = C^H^^kl + klt) + C 2 {k)H^ /VW ^kl + klt) (8.14) 

Using H\l\z) = — ^(J^e™ - J_ lv ( z )) and H%\z) = ^(-J iu ( z ) e — + J_ iu (z)), 
we can rewrite the general solution as, 



fc (t) = A(k)J ikx/VW (Jk 2 y + kit) + B{k)J_ ikx/VW Ukl + (8.15) 
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Therefore, the WKB solution corresponds to imposing the initial condition A(k) = and 
choosing B(k) appropriately, as J_ ik /^y{^j k^ + k\t) has the same time-dependence at early- 
times as our WKB solution. 



We will be interested in the large momentum regime of non-planar wavenumbers, i.e. 
ki 3> H, implying that the factor (1 — r 2 x(r)) in the denominator doesn't vanish anywhere 
(since both r and x are fractions). In the regime x ~ 1 (i.e. early times), the WKB condition 
obviously holds for any momentum (unless r ~ 1 which is the planar case discussed in 
Appendix B). In [9], it was shown that the condition holds for high- momentum modes as long 
as jj 3> expa(i). Therefore, the time at which the WKB solution should be matched with 
the late-time de Sitter solution has a natural choice, r*, such that, 

e aM = ^ (8-16) 

The above equation implicitly states that r* corresponds to late times when e a ( T *) 3> 1. In 
terms of real time, this condition implies that e Q (**) « e #/** ^> i ) Q r t* > U so ~ -gj. Also, if 
t e denotes the time of horizon exit for a given mode of wavenumber k, we have e Hlte = jj. 
This suggests a simple relation between t* and t e , viz. 



t e « 2U (8.17) 

Returning to the problem of matching the modes - since k 3> Hj, x(t*) ~ exp (— 6a(r*)) ~ 
0. Therefore the WKB solution around r = r* is given by expanding equation (8.7) around 
x = and then plugging in the values of x and the frequencies at t = t* . The solution can 
be expanded in powers of e = \J~^~ and one needs to retain terms to the order at which the 
direction dependence first appears. It turns out that it is sufficient to retain terms up to the 
order e 3 and to this order x(t*) and the frequencies are given as, 

x(Q = 2{H I /kf/\l - (Hj/kf/ 2 ) (8.18) 

- 2 2\r H I^ 



k 2 
Hi 



1 + ( §" r2)( T )3/2 



(8.19) 



= w 2 L _ 2 #A (8 20) 

Now, in the de Sitter regime the solution to equation (7.5) is given as, 

<j) k = A+Mr) + ^4-0- (r) (8.21) 
where the modes 4>± are given as, 

Mr) = (±i + ^{-^rfl^j exp (±^-(-3^^ (8.22) 
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Matching the de Sitter solution with the WKB solution at r = r*, we obtain the following 
equations for A + and A-, 



A+<// + (t*) + A_0'_(t*) = 



A + </>+(t*) + A_0_(r*) 



1 




^ f 1 - ^ 
" s k \ 2k 



i 1 + 



2A; 



1 

2w* 



(8.23) 



(8.24) 



where we have absorbed an overall phase in the definition of A + and A_. Solving for A± 
from the above equations, we have, 



A+ = 



2^2H} 



(2-e 2 )+2ier--l)+0(e*) 



exp 



A- = - 



2 v / 2#7 LV 3 



_ r 2 e 3 + 0(£ 4 



exp 



One can easily verify that these coefficients obey the normalization condition, 



(8.25) 
(8.26) 

(8.27) 



In the main body of the paper, we have redefined A± — 
condition reads, 

\A,\ 2 - \A\ 2 = 1 



H 



2^5 A±, such that the normalization 

(8.28) 



9. Appendix B : Planar Modes 

In this section, we discuss the solution of the classical equation of motion for the planar 
modes of a scalar field evolving in the axially symmetric Bianchi I background specified in 
Appendix A. The computation is analogous to the non-planar case and hence we only provide 
a summary of the main results. For details, the reader may refer to [10]. 
As before, we start from 

+u J (r) 2 ) ( f>k = (9.1) 



dr 2 



where the frequency squared is given as oj(t) 2 = 2 ^^f (1 — r 2 x) with r 2 = kv t^ z and x(r) = 

a; 3 

1 — e 6HlT = exp (— 6a)(y / exp (6a) + 1 — 1) and x(t) varies from 1 to as r changes from — oo 
to 0. 

Now, recall that the WKB parameter is given as, 



e = 



1 - x(t) 



k (x/2)V3(i _ r 2a;( r ))i/2 ^l-r 2 ^) 



+ 1 



(9.2) 



For k x <C k y ,k z , the denominator (1 — r 2 x(r)) 1 / 2 can be an extremely small fraction 
at early times (when x(r) — > 1) and the WKB parameter exceeds unity. Therefore, in this 
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regime one needs to solve the equation above directly and impose some initial condition on 
the modes. This is subsequently matched at later times with the WKB solution, which, in 
turn, is matched to the late time de Sitter solution. 

At early times, the frequency squared can be decomposed as to 2 ~ 23 (k 2 + 5uj 2 e® HlT ), where 
5uj 2 « k 2 + k 2 . the equation of motion in the limit r — > — oo therefore has the general solution, 

(ft = Ci(fc) J- iakx/Hl (qe 3HlT ) + C 2 {k)J iakx/Hl {qe m n (9.3) 



with q = a\/5uj 2 /Hi = aJk 2 + k 2 /Hj ( a = is a numerical constant). 



y 1 v 2/ -* v 3 
Following [10], we fix the initial condition as, 



In the limit r — > — oo, this solution assumes the form of an incoming waveform. By noting 
that T(l + iv) = N /3^, we have 



* T -* -»» - feS^ r(i - ik/W ^' 72 '"^' (9 - 5) 

= — == exp [—3iak x T + iip] 
y6ak x 

where ip is a time-independent phase factor defined as e 1 ^ = {a^jk 2 + k 2 /Hj)~ mkx / Hl . 

At later times, the WKB condition becomes valid again and the solution to the equation 
of motion is given as, 



" - : exp [— i / dr cD] H — ^= exp [i I dr Cj] (9-6) 



2Cj J Tn v 2w 



ro V */ ro 



where tl> has to be determined from the equation Gj 2 = u 2 — \ {^j- — ) ; as before. 
Let ro be the time at which and 0| are matched to determine -B-t and this has to be chosen 
such that both the early time solution (equation 9.4) and the WKB solution are valid. A 
suitable choice is given by, 



o 3Hjt - 1 / H i .. 1 I Hi 



«Y {5u 2 )2 a\ k 



The matching procedure leads to the following solution for B±, 



B+n — e v/ V H i ' (9.8) 



B_ « =e v ' V */ ' (9.9) 

y^l _ e -2nak x /H^ 
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Finally, we match the WKB solution to the de Sitter solution 



+ im A - {-< + |(- 3ff ' T > 1/3 ) exp ("I <- 3ff ' T > 1/3 ) < 9 - 10 > 



at a matching time r = r*, such that e a ( r *) = y 777 and this leads to the following solution 
for A + and A-, 



A + = 



yj\ _ e -2n(ask/Hi) 



f ,-n{ask/H I ).-i(n/4:-bk/H I ) 

A_ = - — (9.11) 

J\ _ e -27r(asfc/.ffj) 

where s = k x /k <C 1 is a measure of "planarity" of the mode while a = and b = 
2 are numerical constants. 

The above expressions for A± completely specify the late-time vacuum state of the 
theory. 2-point and the 3-point correlation functions of the scalar field in the "in-in" formal- 
ism can now be computed with respect to this vacuum. 
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